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Abstract 



We consider the correlation structure of the random coelTicients for a 
class of wavelet systems on the sphere (labelled Mexican needlets) which 
I was recently introduced in the literature by [15] . We provide necessary and 

sufficient conditions for these coefficients to be asymptotic uncorrelated 
in the real and in the frequency domain. Here, the asymptotic theory is 
_ developed in the high frequency sense. Statistical applications are also 

^ . discussed, in particular with reference to the analysis of Cosmological 

' data. 
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1 Introduction 

There is currently a rapidly growing literature on the construction of wavelets 
systems on the sphere, see for instance [H], [2], [7], [12], [5D|, [3S] and the 
references therein. Some of these attempts have been explicitly motivated by 
extremely interesting applications, for instance in the framework of Astronomy 
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remarks. The research of Xiaohong Lan was supported by the Equal Opportunity Committee 
of the University of Rome Tor Vergata. 
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and/or Cosmology; concerning the latter, special emphasis has been devoted to 
wavelet techniques for the statistical study of the Cosmic Microwave Background 
(CMB) radiation ([15]). CMB can be viewed as providing observations on the 
Universe in the immediate adjacency of the Big Bang, and as such it has been 
the object of immense theoretical and applied interest over the last decade [5]. 

Among spherical wavelets, particular attention has been devoted to so- 
called needlets, which were introduced into the Functional Analysis literature 
by [551 mi ; their statistical properties were first considered by [5] , [H , [13] and 
[ITj . In particular, it has been shown in [31 that random needlet coefficients 
enjoy a capital uncorrelation property: namely, for any fixed angular distance, 
random needlets coefficients are asymptotically uncorrelated as the frequency 
parameter grows larger and larger. The meaning of this uncorrelation property 
must be carefully understood, given the specific setting of statistical inference 
in Cosmology. Indeed, CMB can be viewed as a single realization of an isotropic 
random field on a sphere of a finite radius ([9]). The asymptotic theory is then 
entertained in the high frequency sense, i.e. it is considered that observations 
at higher and higher frequencies (smaller and smaller scales) become available 
with the growing sophistication of CMB satellite experiments. Of course, uncor- 
relation entails independence in the Gaussian case: as a consequence, from the 
above-mentioned property it follows that an increasing array of asymptotically 
i.i.d. coefficients can be derived out of a single realization of a spherical random 
field, making thus possible the introduction of a variety of statistical procedures 
for testing non-Gaussianity, estimating the angular power spectrum, testing for 
asymmetries, implementing bootstrap techniques, testing for cross-correlation 
among CMB and Large Scale Structure data, and many others, see for instance 
[3], [3, [HI, [H], [13], [29], [23], [8], [5], [20], [3l]. We note that the relevance 
of high frequency asymptotics is not specific to Cosmology (cf. e.g. financial 
data). 

Given such a widespread array of techniques which are made feasible by 
means of the uncorrelation property, it is natural to investigate to what extent 
this property should be considered unique for the construction in [151 HZ] j or 
else whether it is actually shared by other proposals. In particular, we shall 
focus here on the approach which has been very recently advocated by [TS] , see 
also [14] for a related setting. This approach (which we shall discuss in Section 
2) can be labelled Mexican needlets, for reasons to be made clear later. Its 
analysis is made particularly interesting by the fact that, as we shall discuss 
below, the Mexican needlets can be considered asymptotically equivalent to the 
Spherical Mexican Hat Wavelet (SMHW), which is currently the most popular 
wavelet procedure in the Cosmological literature (see again [25l). As such, the 
investigation of their properties will fill a theoretical gap which is certainly of 
interest for CMB data analysis. 

Our aim in this paper is then to investigate the correlation properties of 
the Mexican needlets coefficients. The stochastic properties of wavelets have 
been very extensively studied in the mathematical statistics literature, starting 
from the seminal papers [TOl[Tll. We must stress, however, that our framework 
here is very different: indeed, as explained earlier we shall be concerned with 
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circumstances where observations are made at higher and higher frequencies on 
a single reahzation of a spherical random field. As such, no form of mixing 
or related properties can be assumed on the data and the proofs will rely more 
directly on harmonic methods, rather than on standard probabilistic arguments. 

We shall provide both a positive and a negative result: namely, we will 
provide necessary and sufficient conditions for the Mexican needlets coefficients 
to be uncorrelated, depending on the behaviour of the angular power spectrum 
of the underlying (mean square continuous and isotropic) random fields. In 
particular, on the contrary of what happens for the needlets in [551 HZ] i we shall 
show that there is indeed correlation of the random coefficients when the angular 
power spectrum is decaying faster than a certain limit. However, higher order 
versions (already considered in [15]) of the Mexican needlets can indeed provide 
uncorrelated coefficients, depending on a parameter which is related to the decay 
of the angular power spectrum. In some sense, a heuristic rationale under 
these results can be explained as follows: the correlation among coefficients is 
introduced basically by the presence in each of these terms of random elements 
which are fixed (with respect to growing frequencies) in a given realization of 
the random field, because they depend only on very large scale behaviour (this 
is known in the Physical literature as a Cosmic Variance effect). Because of 
the compact support in frequency in the needlets as developed by [^S] [17] , these 
low-frequency components are always dropped and uncorrelation is ensured. On 
the other hand, the same components can be dominant for Mexican needlets, in 
which cases it becomes necessary to introduce suitably modified versions which 
are better localized in the frequency domain (i.e., they allow less weight on very 
low frequency components). 

As well-known, there is usually a trade-off between localization properties 
in the frequency and real domains, as a consequence of the Uncertainty Prin- 
ciple ("It is impossible for a non-zero function and its Fourier transform to be 
simultaneosuly very small", see for instance [18j). In view of this, an interesting 
consequence of our results can be loosely suggested as follows: the better the 
localization in real domain, the worst the correlation properties. This is clearly 
a paradox, and we do not try to formulate it more rigorously from the mathe- 
matical point of view - some numerical evidence will be collected in an ongoing, 
more applied work. However, we hope that the previous discussion will help 
to shed some light within the class of needlets; in particular, it should clarify 
that the uncorrelation property of wavelets coefficients does not follow at all 
by their localization properties in real domain. Indeed, given the fixed-domain 
asymptotics we are considering, perfect localization in real space does not en- 
sure any form of uncorrelation (all random values at different locations on the 
sphere have in general a non-zero correlation). 

The plan of the paper is as follows: in Section 2 we shall review some basic 
results on isotropic random fields on the sphere and the (Mexican and stan- 
dard) needlets constructions. In Section 3 and 4 we establish our main results, 
providing necessary and sufficient conditions for the uncorrelation properties to 
hold; in Section 5 we review some statistical applications. 
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2 Isotropic Random Fields and Spherical Needlets 
2.1 Spherical Random Fields 

In this paper, we shall always be concerned with zero-mean, finite variance and 
isotropic random fields on the sphere, for which the following spectral represen- 
tation holds, in the mean square sense: 

T{x) = ^ airnYim{x) , X S S"^, (1) 
Irn 

where {0;™}/^, = —I,..., I is a triangular array of zero-mean, orthogonal, 
complex- valued (for to ^ 0) random variables with variance E|a;mP — Gi, the 
angular power spectrum of the random field. The functions {Yim{x)} are the 
so-called spherical harmonics, i.e. the eigenvectors of the Laplacian operator on 
the sphere [H], [M], [32 



sm if 



-l{l + l)Yi^i^,^) 



Ylrn{i},ip) 



where we have moved to spherical coordinates x — {'d,(p),0<-d<TT and 
< if < 2tt . It is a well-known result that the spherical harmonics pro- 
vide a complete orthonormal systems for There are many routes for 
establishing ([1]), usually by means of Karhunen-Loeve arguments, the Spectral 
Representation Theorem, or the Stochastic Peter- Weyl theorem, see for instance 
[T]. The spherical harmonic coefficients aim can be recovered by means of the 
Fourier inversion formula 



aim ^ / T{x)Yim{x)dx . (2) 



If the random field is mean-square continuous, the angular power spectrum {C'l} 
must satisfy the summability condition 

^{21 + l)Ci < 00 . 
I 

We shall introduce a slightly stronger condition, as follows (see [51 HI fT51 121) ). 

Condition 1 For all B > 1, there exist a > 2, and {.9j(-)}j=i 2 sequence 
of functions such that 

Ci = l-"9A^) > for B^~' < I < B^+i , j = 1, 2, ... (3) 



wher 



Cq < gj < Co for all j £ N , and sup sup |- — gj{u)\ < Cr , 

3 B-^<u<B du'' 

some Co, ci, ...cm > , 7\f e N . 
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In practice, random fields such as CMB are not fully observed, i.e. there 
are some missing observations in some regions of S*^; ([2]) is thus unfeasible in 
its exact form, and this motivates the introduction of spherical wavelets such as 
needlets. 



2.2 NPW Needlets 

The construction of the standard needlet system is detailed in [551 HZ] ! we can 
label this system NPW needlets and we sketch here a few details for complete- 
ness. Let (p he a C°° function supported in |^| < 1, such that < <j){^) < 1 and 
= 1 if < l/B, B>1. Define 

b'iO = '^(|) - m > so that viei > 1 , E Jj) = 1 • (4) 

It is immediate to verify that 6(^) ^ only if -g- < |C| < B. The needlets frame 
\^ipji^{x)} is then constructed as 

/ ' 

fjki^) V^E^(;^) E Yirniijk)YiM ■ (5) 

I m— — l 

Here, {Xjk} is a set of cuhature weights corresponding to the cuhature points 
{£,jk} ; they are such to ensure that, for all polynomials Qi{x) of degree smaller 
than B^+i 



^Qi{^jk)^3k^ I Qi{x)dx . 
J 



The main localization property of is established in 26j, where it is 

shown that for any M e N there exists a constant cm > s.t., for every ^ £ 

More explicitly, needlets are almost exponentially localized around any cuha- 
ture point, which motivates their name. In the stochastic case, the (random) 
spherical needlet coefficients are then defined as 



j^ T{x)ipjk{x)dx = y%^^b{^) E a;™i^m(Cjfe) ■ (6) 



I m— — l 

It is then immediate to derive the correlation coefficient 



i>i 
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where Pi is the ultraspherical (or Legendre) polynomial of order i and degree 
I; the last step follows from the well-known identity [34| 

m— — l 

The capital stochastic property for random needlet coefficients is provided by 
[5 , where it is shown that under Condition [1] the following inequality holds 

C 

\Corr (/3 fe,/3 < — g-, some Cm > , (7) 

where d (CjTcCjfe') = arccos ((Cjfej ^jTc')) is the standard distance on the sphere. 



2.3 Mexican Needlets 

The construction in [15] is in a sense similar to NPW needlets in [26l [21] (see 
also [l^), insofar as a combination of Legendre polynomials with a smooth 
function is proposed; the main difference is that for NPW needlets the kernel is 
taken to be compactly supported, which allows on one hand for an exact recon- 
struction function (needlets make up a tight frame), at the same time granting 
exact localization in the frequency-domain. It should be added, however, that 
the approach by ^15j enjoys some undeniable strong points: firstly, it covers 
general oriented manifolds and not simply the sphere; moreover it yields Gaus- 
sian localization properties in the real domain. A further nice benefit is that it 
can be formulated in terms of an explicit recipe in real space, a feature which 
is certainly valuable for practitioners. In particular, as we report below in the 
high-frequency limit the Mexican needlets are asymptotically close to the Spher- 
ical Mexican Hat Wavelets, which have been exploited in several Cosmological 
papers but still lack a sound stochastic investigation. 

More precisely, propose to replace b{l/B^) in © by f{l{l+l)/B'^^), where 
/(.) is some Schwartz function vanishing at zero (not necessarily of bounded 
support) and the sequence + 1)};=! 2 represents the eigenvalues of the 

Laplacian operator A52. In particular, Mexican needlets can be obtained by 
taking /(s) = sexp(— s), so to obtain 

The resulting functions make up a frame which is not tight, but very close to, 
in a sense which is made rigorous in [15 . Exact cubature formulae cannot hold 
(in particular, {Xjk} are not exactly cubature weights in this case), because 
polynomials of infinitely large order are involved in the construction, but again 
this entails very minor approximations in practical terms. More generally, it 
is possible to consider higher order Mexican needlets by focussing on /(s) = 
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s^exp(— s), so to obtain 

i>i 

The random spherical Mexican needlet coefficients are immediately seen to be 
given by 

l^jk:,p = / T{x) V'jfc {x)dx^ j X! X! (2^) V'ife;p (a;) dx 

i>l m=-l 

whence their covariance is 

Throughout this paper, we shall only consider weight functions of the form 
/(s) = sPexp(— s). It is certainly possible to consider more general construc- 
tions; however this specific shape lends itself to very neat results, allowing us to 
produce both upper and lower bounds for the coefficients' correlation. Also, it 
makes possible a clear interpretation of the final results, i.e. the effect of varying 
p on the structure of dependence is immediately understood; this is, we believe, 
a valuable asset for practitioners. In the sequel, we shall drop the subscript p, 
whenever possible without risk of confusion. 



Remark 2 As mentioned earlier, it is suggested from results in fJlSf ) that Mex- 
ican needlets provide asymptotically a very good approximation to the widely 
popular Spherical Mexican Hat Wavelets (SMHW), which have been used in 
many physical papers; the asymptotic analysis of the stochastic properties of 
SMHW coefficients is still completely open for research. The discretized form of 
the SMHW can he written as 

where the coordinates y — 2 tan | follows from the stereographic projection on the 
tangent plane in each point of the sphere; here we take 9 — Ojk {x) :— d{x,^jf.). 
Now write 

by following the arguments in '151 and developing their bounds further, it can 
be argued that 

\^MO;B-^) - K,ki^^,.^^{e)\ = B-W {mm{d''B'^ ,1}) , (8) 
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for some suitable normalization constant Kjk > 0. Equation (0) suggests that 
our results below can be used as a guidance for the asymptotic theory of random 
SMHW coefficients. The validity of this approximation over relevant cosmolog- 
ical models and its implications for statistical procedures of CMB data analysis 
are currently being investigated. 

3 Stochastic properties of Mexican needlet co- 
efficients, I: upper bounds 

As mentioned in the Introduction, having estabhshed ([7]) opened the way to 
several developments for the statistical analysis of spherical random fields. It is 
therefore a very important question to establish under what circumstances these 
results can be extended to other constructions, such as Mexican needlets. In 
this and the following Section, we provide a full characterization with positive 
and negative results. We start by writing the expression for the correlation 
coefficients, which is given by 

E(^F(^Fe-'('+i)(^-^"+^-^^^)(2/ + l)GP,((C,,,„C,..J) 
1>1 

I E (^)4Pe-2'('+i)/s^^'^ [21 + l)Cz I I S (^)4pg-2;(i+i)/s^J2 + i)q | 

We now provide upper bounds on the correlation of random coefficients, as 
foUowfl 

Theorem 3 Assume Condition Q] holds with a < 4p + 2 and M > Ap + 2 — a; 
then there exist some constant Cm > 



\Corr {Pjikiipj (^j2k2;p) 



< 



Cm 



(1 + B(n+T2)/2-ios,(n+n)/2d ' 

(9) 



Proof. We prove ([9]) following some ideas in [26] . For notational simplicity, we 
focus first on the case where ji = J2; we have 

Corr (/3,i^.„,/3,fe,.J = -— ■ . (10) 

V 3k,pi i^jL ,pj E ( 2I )2Pc~2'^^-B^"' 2i+l Q ^ ' 

1>1 ^ 



^ While finishing this paper, we learned by personal communication that working indepen- 
dently and at the same time as us, A.Mayeli has obtained a result similar to Theorem |3] for 
the case J 1 = j2 1 see |24| . The statements and the assumptions in the two approaches are not 
equivalent and the methods of proofs are entirely different; we believe both are of independent 
interest. 
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Now replace Ci = I "'gj (57) in the denominator of the above representation, 
for which we get 

Denoting 9 = arccos k^^j k')i the numerator can be written as 

^ A. "^'ttX fcos^-cos^)^/^ ^ ^ 



where we have used the Dirichlct-Mehler integral representation for the Legendre 
polynomials [16j . 

The following steps and notations are very close to [26] . We write 



' 47r SJ / V 2 , 

i>i \ / \ / 

where 

By Poisson summation formula, we have 

l>l /GZ /jGZ 

Denote 

G„^, (0 e^-^g, it) e-2*(*+^-^)/(5_,_^) . (13) 
Let us now recall the following standard property of Fourier transforms: 

Some simple computations yield 



47r 



+ ^ B-i^'^'^^ + 2^-^ \ (.) |.=B.(2..^.) , 



9 



where 

For all positive integers k < M, we can obtain 



oo ,k 



< 



fcr(m+2p-a)c, ^_,(2p+i+„.-o-fc) for 2p + 1 - a + r 
kV {m + 2p- a) Cg{j log B) , for 2p + 1 - a + m 

(14) 



where Cg = max{co, cm} and 



L{p, m, a, k) 



r {2p + I + m ~ a ~ k) when (2p + 1 + to - a - fc) > , 
(r (A; - 2p - 1 - TO + a))~^ when {2p + 1 + m - a - k) < 



It should be noticed that our argument here differs from the one in |26j , because 
wc cannot assume the integrand function on the left-hand side to be in for 
all k < M. Let us now focus on the case where fc = 4p + 2 — a, with m = 2p+l; 
we obtain 



< 



oo ik 



\B^ (27r^ - if) 



4p+2-a 



therefore 



where 



< 



< (4p + 2 - a) r (4j5 + 1 - a) Cg(j log B); 



(27r/i - (^) 



4p+2-a 



C2p+i,a.g - (4p + 2 - a) r (4p + 1 ~ a) log S . 
The modifications needed for other cases are obvious and we obtain 



_g-(2p+l-m)j 



< 



G R-i(4p+2-a) 
{2TTfl - Lp) 



Ap+2-a 



where 



_ (4p + 2 - a) r (to + 2p - a) Gg 
L{p,'m,a,k) 

Now let Ga^g — Taax{Gm,a,g, m — 0, 2p + 1} ; we have 
h± (27r/i) 



< 



< 



5(2-a)i 

An 



2p 



G, 



{in=Q \ ^ 



2p\ i3-j(4p+2-Q) j^-j{ip+2-a)+logg j 



2'^PGa,gB-^P3+^°SB 3 

{2Trn T <y5) 



ip+2-a ' A* = 1; 2, ... 



(15) 
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Therefore 

Ir'r, I < o'^PC I ^ L \^ ]. 1 R-4pj+iog^i 

^ 2^" ( + |4p + 1 - a| 7r"-4f-i) 

Hence, for the numerator of the correlation we have the bound 

_ 4 -^i ■ r (2^^ + |4p + 1 - a| ^"-4p- A 

(HI]) < CS-4"+'°sbJ / J^Jf Ld^ , 

Je {cos9 — cos(f) ' 

As in [Mj, when O<0<7r/2, we can get the following inequality 

If 7r/2 < < TT, letting 9 = n — 9,ip = n — ip, we can obtain the same bound. 
Going back to , we obtain 



Corr{(i^,,.^(i^^,,) < C 

< (^g,a-4p-2^-0-logB j)(4p+2-a) ^ aS j ^ OO , 

We thus get inequality © for j'l = j2 ■ 

Now let us consider ji 7^ As the proof is very similar to the arguments 
above, we omit many details. In the sequel, for any two sequences ai,bi, we write 
o; « bi if and only if a; = 0{bi) and bi = 0{ai). 

First, we consider the variance of the random coefficients, which can be 
represented by: 



^(fcil )4j>e-2/(i+i)/B^^i (2Z+1)G « I (^y' + j^^ t^P+^- 



Qg-2t(t+B-Ji) 



9ji (t) dt 



\ 4p + 2 — a ' ' 

therefore 

1/2 r ^ 1/ 



J>1 J [i>l 



^ 5(l-"/2)0l+j2) J ^-ii(4p+2-a) , (5 Q\ I J '^0 g72(4p+2-a) i Q (1) I 

\4p + 2-Q; l4p+2-a 



11 



Without loss of generality, we can always assume ji < j2- We can implement 
the same argument as before, provided we replace CB,gj in (fT2l) by 



WjYl ('^nj2 + {l) - (0) 



2i' 
i>i 



where 



Again, by Poisson summation formula, we have 



i>i 
Denote 

by the same argument and notation as in (|15p . we have 



47r 



m— 1 

^2p(ji-j2) + (2-Q)iig±iv r rf2p+l 



47r [ (iw2p^ 

Therefore 



It is then straightforward to conclude as in the case where ji ^ j2 , to obtain 



^Q-4p-2 ^-2p0i+j2)+logB jl 
^ (jQa-'ip-2^-[(ji+j2)/2-logsUl+j2)/2]{4p+2~a) ^ aS J2 ^ OO 



Thus dSl) is established. 
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Remark 4 By careful manipulation, we can obtain an explicit expression for 
the constant Cm in (0^, i-e. 

Cm = 22P^*^+1 (4p + 2 - af T (4p + 1 - a) coCg log B. 

The previous result shows that Mexican needlets can enjoy the same uncor- 
relation properties as standard needlets, in the circumstances where the angular 
power spectrum is decaying "slowly enough" . The extra log term in ([5]) is a con- 
sequence of a standard technical difficulty when dealing with a boundary case 
in the integral in 



Remark 5 To obtain central limit results for finite- dimensional statistics based 
on nonlinear transformations of the Mexican needlets coefficients, it would he 
sufficient to consider the case where j = f . The asymptotic uncorrelation we 
established in Theorem\^is stronger; indeed, for many applications it is useful to 
consider different scales {j} at the same time. Because of this, it is also impor- 
tant to focus on the correlation of Mexican needlet coefficients at different j, j' . 
We stress that the need for such analysis was much more limited for standard 
needlets; indeed in the latter case, given the compactly supported kernel 5(.), 
the frequency support of the various coefficients is automatically disjoint when 
\j — j'\ > 2, whence (for completely observed random fields) standard needlet 
coefficients can be correlated only for \j — j'\ = 1. 

4 Stochastic properties of Mexican needlet co- 
efficients, II: lower bounds 

In this Section, we complete the previous analysis, establishing indeed that the 
random Mexican needlets coefficients are necessarily correlated at some angular 
distance in the presence of faster memory decay. This is clearly different from 
needlets, which are always uncorrelated. As mentioned in the Introduction, the 
heuristic rationale behind this duality can be explained as follows: it should be 
stressed that we are focussing on high-resolution asymptotics, i.e. the asymp- 
totic behaviour of random coefficients at smaller and smaller scales in the same 
random realization. For such asymptotics, a crucial role can be played by terms 
which remain constant across different scales. In the case of usual needlets, 
which have bounded support over the multipoles, terms like these are simply 
dropped by construction. This is not so for Mexican needlets, which in any 
case include components at the lowest scales. These components are dominant 
when the angular power spectrum decays fast, and as such they prevent the 
possibility of asymptotic uncorrelation. In particular, we have correlation when 
the angular power spectrum is such that a > -\- 2. 



13 



Theorem 6 Under condition{J\ for a > 4p + 2, V e G (0,1), there exists a 
positive < e [1 + Cq j sucri triat 



/il^ inf Corr {Pjk-p, Pjk'-p) > 1 - e ; (16) 
for all {ijk^ijk'] such that rf(Cj7c,Cjfe') < ^ ■ 

Proof. We first divide the variance of the coefficients into three parts, as follows 

E + E + E ) (^)'%-<->"'"?i±ic, 



Intuitively, our idea is to show that the first sum is of exact order 0{B^^ x 
^(a-4p-2)j) ^ 0{B-^Pi), while the second two are smaher (0(B(2-")j) and 
o{B^^^°'^^), respectively). Indeed, for the first part we obtain easily 

l</<eiSi ^ ^ l<;<ciBJ 

< 2^B(2-"b / 2;4p+i-"d:r = 2^ ^ ^ ^(s-ab^ 

TT 7r(a-4p-2) 

and 

K^tr^sA / 47r 2^co(a-4p-2) 

Similarly, for the second part 



( 



4p+2-a 4p+2-Q , „ 



-2e^n(2-a)j < ( + '^) \ -2l{l+l) / B^^ '^[±±^ 

TT (a — 4p — 2) 

and for the third part 

2p 



\ B^3 J At: Att i„ 
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The last inequality follows from the asymptotic formula 



Jy ^ y 



We have then that 

^ _ J /[([+l)V%_2,(i+l)/B=.-2[+l^' 
I l<l<eiB3 



1 1 ji +^) g-2i(i+i)/B2^ 2/ + 1 ^^ 



l>f.2B3 



= 0(Bi(4p+2-")) ^ p(^) , as i ^ 00 . 

On the other hand, for any positive e < 1, if we choose ei = NB~^ , where A'' is 
sufficiently large that 7V4p+2-a (i + 1^) < 1, we obtain that 



Co £ 

whence 

< - as 7 ^ oo . 
2 

Thus we obtain that 



> _ (et^'+'-« - e^^+'-«) CO 



2i , 



-1 



> cB^^P\ some c> . 

More explicitly, the variance at the denominator has the same order as for the 
summation restricted to the elements in the range \ <l < t\BK 
To analyze the numerator, we start by recalling that 

sup Pi{cos9) = P;(cosO) = 1 , and sup |-^Pi(cos6')| < 3/ . 

6ie[0,7r] 6>e[0,7r] "'(' 

As a consequence, for any £ > there exists a. 6 > 0, s.t. if < ^ < 5 < e/6N, 
then, 

|P;(cos6') - P( (cos 0)1 < 3/6' < e , 
for all I > N, where the above inequalities follow from 

n 

0<cos(6'o + 6')-cos6'o = 2sin2- < 61 , V 6»o € [0,27r) . 
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Therefore, for any £,jk,S.jk' G ^"^ ^.t. arccos (f^^,, i^^j,,) < S, we have 



Corr (/3jfc,/3jfe, 



1>1 



2p 

(,-21(1 + 1) /B^i '21+1 (J 



> 



2p 

g-2((i + l)/S2j 2i+l(7 



E(^) 



1>N 



e(^)" 



+o(l) = (l-e') + o(l) , asj->oo , e'- 



l + e/2 ^ ' ^ ^ ^ ^ ' 2 + £ 

Thus the proof is completed. ■ 

As mentioned ear her, the results in the previous two Theorems illustrate 
an interesting trade-off between the localization and correlation properties of 
spherical needlets. In particular, we can always achieve uncorrelation by choos- 
ing p > (a — 2)/4; of course a is generally unknown and must be estimated 
from the data (in this sense standard needlets have better robustness proper- 
ties). Introducing higher order terms implies lowering the weight of the lowest 
multipoles, i.e. improving the localization properties in frequency space. On the 
other hand, it may be expected that such an improvement of the localization 
properties in the frequency domain will lead to a worsening of the localization in 
pixel space, as a consequence of the Uncertainty Principle we mentioned in the 
Introduction (see for instance [H]). We do not investigate this issue here, but 
we shall provide some numerical evidence on this phenomenon in an ongoing 
work. 
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Remark 7 In Theorem\^ we decided to keep the assumptions as close as pos- 
sible to Theorem O in order to ease comparisons and highlight the symmetry 
between the two results. However, it is simple to show that the correlation result 
holds in much greater generality, for angular power spectra that have a decay 
which is faster than polynomial. In particular, assume that 

Ci = H{l)exp(-F) ,1 = 1,2,... 

where H(l) is any kind of polynomial such that Hil) > c > and p > 0. Then it 
is simple to establish the same result as in Theorem O by means of a simplified 
version of the same argument. The underlying rationale should be easy to get: 
for exponentially decaying power spectra the dominating components are at the 
lowest frequencies, and they introduce correlations among all random coefficients 
which cannot be neglected. 

5 Statistical Applications 

The previous results lend themselves to several applications for the statistical 
analysis of spherical random fields, in particular with a view to CMB data 
analysis. Similarly to [3], let us consider polynomials functions of the normalized 
Mexican needlets coefficients, as follows 

^ Q 

where IIq{.) denotes the g-th order Hermite polynomials (see [33]), Nj is the 
cardinality of coefficients corresponding to frequency j (where we take {^jfe} to 
form a i?~^-mesh, see so that Nj w B^^), and {wuq} is a set of deterministic 
weights that must ensure these statistics are asymptotically nondegenerate, i.e. 



Condition 8 There exist jo such that for all j > jo 

rank{VLj) = U ,VLj := Eh^^h'^. , h^^ := {hi^N,, ■■■,hu,N,)' ■ 

Condition ([8]) is a standard invertibility assumption which will ensure our 
statistics are asymptotically nondegenerate (for instance, it rules out multi- 
collinearity) . Several examples of relevant polynomials are given in [3]; for 
instance, given a theoretical model for the angular power spectrum {C;} , it is 



jk;p 



, u= l,...,U 



E0 



jk;p 
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suggested in that reference that a goodness-of-fit statistic might be based upon 



^^■J fe=l V^^J k=l 



i>i 



The statistic 



r,- = 

iV..- 

fc=i 



can then be viewed as an unbiased estimator for 

i>i 

We refer to [13], [H] for the analysis of this estimator in the presence of miss- 
ing observations and noise, and for its apphcation to CMB data in the stan- 
dard needlets case. Our resuhs below can be viewed as providing consistency 
and asymptotic Gaussianity (for fully observed maps and without noise) in the 
Mexican needlets approach. As always in this framework, consistency has a 
non-standard meaning, as we do not have convergence to a fixed parameter, but 
rather convergence to unity of the ratio Tj/Tj. 

Likewise, tests of Gaussianity could be implemented by focussing on the 
skewness and kurtosis of the wavelets coefficients (see for instance ^25^), i.e. by 
focussing on 

1 r - ^ 1 ^' 



H3if3jk-p) + Hi{f3jk.p)j = —= ^ and 



fc=l y^^3 k=l 



V k=i V i k=i 

The joint distribution for these statistics is provided by the following results: 



Theorem 9 Assume T is a Gaussian mean square continuous and isotropic 
random field; assume also that Conditions [21 [3 are satisfied and choose p > 
{a + i5)/4, some d > 0. Then as Nj oo 

nj^^^hN^ N{o, lu) . 

Proof. The asymptotic behaviour of our polynomial statistics can be estab- 
lished by means of the method of moments. In particular, it is possible to 
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exploit the diagram formula for higher order moments of Hermite polynomial, 
as explained for instance in [28], [33]. The details are the same as in [3], and 
thus they are omitted for brevity's sake. We only note that, in order to be able 
to use Lemma 6 in that reference, we need to ensure that 

C 

\Corr (/?^.,.p,/3^-,.,.p)| < 2+6 ' ™ C> . 

In view of ([9]), this motivates the tighter limit we need to impose on the value 
of p. m 

It may be noted that the covariance matrix flj can itself be consistently 
estimated from the data at any level j, for instance by means of the boot- 
strap/subsampling techniques that are detailed in [3]. Again, the arguments 
of that paper lend themselves to straightforward extensions to the present cir- 
cumstances, as they rely uniquely upon the covariance inequalities for the ran- 
dom wavelet coefficients. Likewise, the previous results may be extended to 
cover statistical functionals over different frequencies, for instance the bispec- 
trum ( |21l \22\ ). A much more challenging issue relates to the relaxation of the 
Gaussianity assumptions, which is still under investigation. 
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